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Abstract 


• A general theory for steady, oscilla- 
'tory or fully unsteady potential cc'^pres- 
'siblc aerodynamics around complex config- 
'urations is presented. .Usirg the finite- 
leloment method to discratire the space 
■problem, one obtains a set cf differential- 
delay equations in time relating the poton- 
jtial to its normal derivative (on the sur- 
face of the body) v:hich is expressed in 
terms of the generalised coordinates of 
the structure. For oscillatory flov/, 
motion consists of sinusoidal oscillations 
around a steady, subsonic or supersonic 
•flow. For fully unsteady flow, the notion 
i.s assumed to consist of constant subsonic 
*or supersonic speed for time t< 0 (steady 
state) and of small perturbations around 
the steady state for time t>0; the solu- 
tion is obtained in Laplace's domain. 

From the potential, the aerodynamic 
(generalized forces are obtained. Therc- 
:fore the final output is the matrix of 
the aerodynamic coefficients, relating 
ithe generalized forces to. the generalized 
jcoordinates, in the form necessary for 
flutter applications. The theory is em- 
bedded in a computer code, SOUSSA (Steady, 
Oscillatory and Unsteady, Subsonic and 
iSuper sonic Acrodynauaics) , which is 
(briefly described. Numerical results 
iare presented for steady and unsteady, 
jsubsonic and supersonic flows and indi- 
cate that the code is not only general, 
Iflexible, and simple to use but also 
]eccurate and fast. 

i' " 

■1. Introduction 

Presented herein is a general formula- 
tion of steady, oscillatory or fully un- 
steady, subsonic and supersonic potential 
‘aerodynamics for an aircraft having arbi- 
trary shape. The objective of this formu- 
ilation is to described the time functional 
jrelationship between aerodynamic potential 
|and its normal derivative (normal wash) in 
a form which can be used for computational 
janalysis. The finite-element method is 
.used for space discretization. The matrix 


!of the aerodynamic influence coefficients, 

'an necuscary foi flutter calculations, is 
[then obtained. Results obtained with the 
(computei program SOUSSA (Steady, Oscilla- 
• tory and Un-Tteady, Suh sonic and Supersonic 
(Aerodynamics) arc also presented. 

i '! 

, The analysis presented heroin is based ; 
'on a new integral formulation, presented in > 
{References 1 and 2, v.’hich includes cor.-.plctc- 
( ly arbitrary motion. However, the ntreri- 
ical implementation (Refs. 3 and 4) was thus 
(far limited to steady and oscillatory flows. 
:0n the other hard, in order to perform a 
! linaar-systcir. analysis of the aircraft, it' 

J is convenient to use more general aerody- 
jna.mic formulations, i.e., fully transient 
(rorponsc for time-domain analysis and the 
j aerodynamic transfer function (Laplar.n 
(transform of the fully unsteady oper.atcr) 

.for frequency-domain analysis. A general 
! formulation for fully unsteady (indlcial) • 
•aerodynamics was presented in Refs. 5 and 6 ' 
•where only very preliminary results vere 
given. Consistent with this type of 
■analysis, the unsteady contribution is as- 
‘sumed to start at time t=0, so that for 
tine t<0 the flow is in steady state. 
Furthermore, consistent with the iinev.r 
flight dynamics analysis, the motion of ( 
the aircraft is assumed to consist of 
small (infinitesimal) perturbations rround 
,the steady-state notion. 

It 'may be noted that within the as- 
sumption of potential aerodynamics, tr.ere 
(exists other methods to evaluate tiio eero- 
. dynamic loads.' Among them, the liftis.g 
surface theories, while flexible and effi- 
cient, are not sufficiently general. On 
'the other hand, finite-element methods, 
though sufficiently general for handling 
complex configurations, are limited to 
steady flows. In addition, they arc usu- 
ally quite cumbersome to use and invari- 
ably require human itervention to define 
; the suitable type of element (source, 

^doublet, etc.) to be used.' For oscilla- 
itory aerodynamics, the doublet-lattice 
method *'* is the only other method, be- 
.sidcs SOUSSA, which can handle subsonic . * 
oscillatory flows around complex configur- 
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ations, while SOUSSA is tlio'only program 
which can atialyzo oscillatory supersonic 
aerodynamicc. 

I Finally, for fully unsteady aorodynam- 
ica, covcral problems have been considered 
'since the initial work by Wagner (P.ef. 10) 
on unsteady inoonprersible two-dincnsionnl 
jflo'./. While rcvi'ral method’s nru available 

I for \;ings in subsonic and supersonic flov/ 
(see Refs. 11 and 12), no other code, be- 
sides SOUSSA, is available for subsonic 
‘and supersonic flows around arbitrary com- 
iplex conflg’urp.tions for either time or 
jfrequency domain analysis. 

The purpose of this paper, is to pre- 
sent recent developments on the formula- 
tion of Ref, 1. In this paper only the 
subsonic formulation is presented in 
details* the supersonic fonrulation is 
only briefly outlined in Appendix A. For 
conciseness, material previously pre- 
sented (in particular, the material of 
Ref. 4) is not repeated herein. 


1 Using the finite-element method to 
I discretize the space problem, one obtains 
•a set of differential-delay equations in 
time relating the potential to its normal 
derivative (on the surface of the body) , 
‘which is expressed in terms of the gener- 
alized coordinates of the structure. For 
oscillatory flow, the motion consists of 
sinusoidal oscillations around a steady, 
subsonic or supersonic flow. For fully 
;unsteady flov;, the motion is assumed to 
consist of constant subsonic or super- 
sonic speed for time t<0 (steady state) 
and of small perturbations around the 
steady state for tine t>0; the solution 
is obtained in Laplace's domain. From 
;the potential, the aerodynamic general- 
ized forces are obtained. Therefore 
the final output is the matrix of the 
aerodynamic coefficients, relating the 
, generalized forces to the generalized 
■coordinates. The theory is embedded in 
a co.-nputer code, SOUSSA (Steady Oscilla- 
tory and Unsteady, Subsonic and Super- 
sonic Aer od ynam ic s) , which is briefly 

[ ''described. Numerical results are presented' 
for steady and unsteady, subsonic and 
supersonic flows and indicate that the 
(code is not only general, flexible, and 
[simple to use but also accurate and fast. 

i 1 


2. Equation for Velocity Potential 


The subsonic aerody.iHm.ic forr.ulation 
used in SOUSSA is briefly presented here. 
IThc supersonic formulation is given in 
I'Appendlx A. Assume the flow to be an in- 
l•finitcsimal perturbation from the steady 
j state flov;. Then standard use of Green's 
I'function method applied to the equation 
•of the velocity potential yields, after 
I linearization, the following inLcgral 
! equation ‘ . 


2n ‘KP„T)=-J^ + 


( 1 ). 



r54» 


3N \ R / 

Ut. 

J R »N 1 B 




a N\ R / 

1^1 

J R BNr’=« 


where is on the surface cf the body, 
VE n ,N is the outer normal 








(2)‘ 


1 X= x/dZ y~y/l Z=z/X T«q^P»A 

I I 

wit h p« J 1 -M and 

i R“[(X-X.)2 + (y_Yj2 + {2„2j2j 


( 3 ) 


’while 


. where 


•I 


li'-TJI.., 


® =M (X-XJ + R 


( 4 ) 


(5) 


I is the time necessary for a disturbance _to 
propagate from P to P*. In addition, 

i is the (open) surface of the wake (k now n 
,i from the steady state solution) and A4» is'^ 
• the potential-discontinuity across the 
■ wake, evaluate_d_ in the direction of the 
'normal, i-. e. if the upper nor- 

mal is us_cd.. It should be noted that the 
!. value of .A']> is not an additional unknown. 


since 


T) » A4-(P^g, T-n) 


( 6 ) 


, where IT is the nondimcnsional time neces-j 
j' sary for the vortex-point to travel (with-i 
- in the steady flow) from the point, i 


TE 


(origin o£~th'e vortex-TTne at tlie trailin’g 
edge) , to the point .P_. For small-pertur- , 
bation steady flov;, II' is given by ‘ I 


n 


B ^ CX - X^g)/M 


( 7 ) 


Equations (1) and (6) fully describe 
• the problem of linearized unsteady sub- 
I sonic potential aerodynamics around 
I complex configurations. In order to ‘solve 
j this problem, it is necessary, in general, 
} to obtain a munerlcal approximation for 
Eq. (1). This is obtained by dividing the 
surface of the aircraft into Ng_ quadri- 
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lateral elements X|, (which nrc described 
in terms of the corner points by use of 
standard finite-ele-ent interpolation 
technique*) and by tssuninq <{' ind to 
be constant within each clement: 


T-e)=='^,^(T- pj,) 


♦ (p, T- 0 )-*h 


( 8 ) 


where 

time dependent values of V and <t> at the 


0 . 


(where 

is the distuzbance-propagation tine 


centroid P^. of S, at the time T - 0. 

n h n 


from to 5j^) . 

Next consider the integrals on the 
wake. In order to facilitate the use of 
Eq. (6 ) , it is convenient to divide the 
wake into strips defined by (steady-state) 
vortex-lines emanating from the nodes on 
the trailing edge. The strips are then [ 
divided into elements Ejj (vr) with nod«=!S j 

along the vortex lines. The potential 
discontinuity is assumed to be cons- 
tant within each element 


A*(P, T -e)“ f'" "®n^ 


(9) 


where (T-9 )• is the value of 

j ^(W) (VV) 

at the centroid P™' of the element e 1-— • 
n , ^ n 

on the wake at time T - 0 (where 0^ is 
the p ropag atio n ti me f rom to . 

Note that according to Eq. (6) 


A* (r)« 


a-n) 

m{n) n 


’ where m>m(n) identifies the traillng- 
pdge point which iscn the same vortex- 

line as the point • Furthermore, 

is the time necessary for the vortex- 
point to be convected from the tralling- 

. edge point to the wake-point P^' . 


It may be worth noting that 


0 , 




where h. and h% 

U 


m 

identify 


"u "A 

the upper and lower trailing-edge nodes 
on the body corresponding to the mth 
node on the trailing-edgo. 

In SOUSSA 9^ - is approximated 
"o "A 

with the value evaluated at the centroids 
of the elements adjacent to the trailing 
edge. This is reasonable in view of the 


* The equation for the eJementn are of 
the type P » Po+C'?i+np2'‘’^n^3 ^ 

< l;-l<r|< 1) . This type of element is 
called hypcrboloidal element and is 
described in details in Ref. 4. 


Kutta condition. ^ ,n > is then evaluated 
from the centroid, as * 




( 11 ). 


with the above approximation, it is 
possible to write * 


A* cry . s .♦ 




h 


nh^h 


( 12 )* 


where * l(S^j^»-i), if h identifioc 
i the upper (lower) point P^ on the body 
j corresponding to the point p^^^ on the 

‘ on 


wake (i.e., near the point 
■ the tnlling edge) , and otJie’r- 

i wise. 


Combining Ecs_^ (1) , (8) and (9) 
I and assuming P»*F| > one obtains 


(13) 










■•■£2 Pr„ s . \ (r -0. -n ) 


n h 


! 

i where 

I 

I 


( 10 ) 


JJt 3 - 


(14) 


P* «= F, 




'P* “ P, 


D.. =-l_ ff J- iL 
2nr JJ^R SN 

^( t )< 


1 . dR - 


|P*“P, 


|P*“Pt 


and 


In “ " ~ ff ~ > 

i 2n JJ-- R 8JM n 

fn 


ih 


P*=P| 


F.“P, 


(15) 


Using the above mentioned hypcrboloidal 
quadrilateral elements the coefficients 
Bjj^, Cjj^, and are evaluated 

analytically with the expressions given 
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• ill* • t'i * 7'< 


ii Ilk »<9f* 


in Ref. 4. The coefficients D.. arc 

jh 

approxinated as D^. =R.. whore 
- jn jn 

"jh- 

Next* taking the Laplace transform with 
zero initial conditions of Eq. (13) yields* 


K] 1^1 "PiJi' 


where and are the Laplace trans- 
forms of**}, ancT*}, -while 




I and 




2 ,k-V‘*°- 




(17) 


nh 


(18) 


• whereas S is the nondimcnsional Laplace 

; parameter.* 


' 3. Boaindary Condition 

i The derivation of the boundary condi- j 
' tion is given in details in Ref. 13. Here 
! the derivation is briefly summarized, 
i The bo'jndary condition obtained by inpos- 
‘ ing that the velocity of the fluid and 
‘ the velocity of the body, v , have the 
; siuce components along the normal to the 
: surface of the body. This yields | 

i 

• • • f 

I ‘t," (v -UooT ) . n (ty'u^ .(19) 


j- where n(t) is the Instantaneous normal 
i to the surface of the body. The unsteady 
I part of the boundary conditions (neglect- 
! ing high order terms) is given by 


v • n + An • T 


( 20 ) 


with An »7 ® 


1 


whore d is the displacem.^t of a point 
of the body. Expressing d v in 

terms of the generalized coordinates 


of q . as 

A 


d 


M 

n**l 


(21) 


* For oscillatory aerodynamics, sotting 

• Aw(T) - and ¥j,(T) « 

yields” the same equation with S«in 


and 




(16) . 


i E % “.ft"* 

n I 

one obtains, in the Laplace domain , 

• - r 

♦■E ■ 

n - l LHs> 

lojXo^lVbs’ ^ 2 ’* / 


( 22 ) 


(23) 


; where » " (23) gives the 

i desired'rcTatiohship between normal wash 
j at any point and tho gonernlized coor- 
■ dinates <j 


j 4. Pressure and C enc -r.nllzed Forces 

[ In order to complete the formulation, 

I the procedure for the evaluation of the 
I aerodynamic pressv.i t' and the generalized 
forces is presented in this section. 

First, consider an averaging scheme 
which imposes that the value of the 
potential at the node p' is the 

average of the val-acs of th^ potential 
! at the centroids of the elements sur- 
j rounding ^ . In other words 


'V 


where is an averaging matrix 

defined as 


(24)1 


8d “ 5 d “ .* 

— • • a, X I - — ,r • 0. XI 

ISjXa^l ^ 8? 


‘kh 


1 




if p; 6 E, 
k h 


(25) 


(i.e., if the P,' is one .of the corner i 

k 


points of the clement £|^) and 


E, . * 0 othervdso 
kn 


(26)' 


In Eq. (25), _is the number of ^ 
I elements which have as one of their . 

corner points. Having evaluated the 
values of at t.ic four corner points 
of each quadrilateral element, the po- 
tential is expressed as 


N' (P) 

h h 


( 27 ). 



ORIGINAL PAGE IS 
OF POOR QUALFIY 


Al . . Cul'i Sll. l* rs 


I I 


I*.. 


I»* •,«' #• I 

' I* |> V* If &I. - 


whcra M' are the firct-ordor ylobal shapo-^ 

functions obtained by asBer.bling local 
shape-functions of the typo ) 


N 


(E): 


1 




(5 +t|^)0l + \) 


(20) 


where “i 1 »hd 1)|^ •+ 1 ere the looa- 
tlon(of~thc cornursi^” of the elec’.cnt 

( C and T)’ are the coordinates over th^ 
eldoent)' so that 


— •ly. — ' 

8S» ‘’dS® 


(29) 


where 2^ 


5 ; 


The pressure coefficient may be evalu- 
ated from the linearized Bernoulli 
theorem as 


VM;jT 9 / 


(30) 


I 


Expressing V^^ in terns of the tangential 
derivatives"^ ♦. and neglecting the con- . 
tribution of the normal cemponent, the j 
pressure coefficient is given by 


- 

B \M hi 


+ I 


'■Sr 


+ 1 


7} 


a* / 


where X- are the contravaniant base 
vectors, with. 3^/32* given by Eq.(29). 

Next consider the generalized aero- 
dynanic forces 




M dE- 
n » 


(32) 


where 


I 


(31) i 


\ . r' 

t- ‘Oj X o'^ • 

« -4(oj xo^2 * p“p. 


(35) 


where ,Oj > o-:d ^2 bhe base vectors of 


the element E,> 
h. 


S. SOOSSA 


where 


The above forni”!) ■ t ion is ir.uJ.' i;'Gntcd 
in the computer prepran SOUESA (Steady, 
Oscillatory and Unsteady, Su'r.-^^onic and . 
Supersonic 7.erodynr::.'ics) . Tr.c-. program 
is an iraprovciucnt of the progran EOSSA 
presented in Ref. 4 and therefore re- 
tains all of the banfe features analyzed 
in Rc£. 4. In par ' ^•ular tl.u program 
besides being goncral and flexible is 
also very simple to use. The only Inputs 
aie the location of the corner points of 
the quadrilateral elements, the Mach number 
’ and the reduced frequency. The wake in 
automatically generated. It should be 
noted that a considerable im;.rrovaincnt • 
with respect to SOSEA'’is that supersonic 
flows are treated exactly the same way 
as the subsonic ones: in particular 
diaphragm.s are not used in SOUS3A4 There- 
fore the basic simplicity in use for the 
subsonic flows is retained i.i supersonic 
flows as well. An additional feature 

which didn't exist in SOSSA is the evalua- 
tion of the aerodynamic influence coeffi- 
cients. Therefore it may bo worth it to 
add a few* comments on the computational 
implementation of the formulation pre- 
sented above. (Details for the following 
expressions are given in Ref. 13.) Note 
first that the relationship between nor- 
mal wash and generalized coordinates, 

‘ (23 ) may be w’rittcn as* 




GO 


(33) 


" " (1) ~ 


(36) 


la the dyneunlc pressure. 

By assuming that the pressure coeffi- 
cient "“e"' is constant within each element 
(consistent with the assumption made on <1> ) > 
Eg. (32) can be expressed as 


• (O-n ) - t Wl.1 • 


(34) 


where 


4' 


j^i 


( 37 ) 


i* In Eqs. (36) to (43) compact matrix 
j notations arc used. Vectors and 
I matrices are underlined. Tildas in- 
I dicate Laplace's transform or eguiva- 
I lent operation. 
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and fj' ' n r.atrix implicitly defined by 
Eq. (23). Fui t'hQi'fr.oro Eq. (16) may bu 
rewritten 




whereas the relationship between pressure 
coefficient and potential may be written 

as 


m 5 


where 


is obtained by combining Eqs. 


(24), (27), (29) and (31). Finally (sea 

Eq. (34)) 


Q -q c (11) 

. - - -P 

or, combining Eqs. (3G) , (39) , (40) and (41), 

Q -qMa 

where the matrix of the aerodynamic inilu-. 
ence coefficients H is given by 

‘ ( 43 , 


Note that and deoend upon 

the nodes 5* ), whereas ^(2) and 

JH' ' do not. The separation of {1 into 
mode-depandent matrices and mode-in- 
dependent matrices is particularly 
useful for automated structural design 
in which the sane geor.etry but different 
nodes are used in each iteration. In 

addition is the only matrix v/hich 

depends in a complicated vray upon the 
complex frequency S. However once the 
coefficients of Eqs. (14) and (15) arc 
evaluated, the evaluation of m(2) re- 
quires only the combination or this 
coefficients according to Eqs. (L7) an^ 

0.8) and the inversion of the matrix Y 
(see eq. 39) . This is particularly 
useful for the evaluation of the aero- 
dynamic influor.ee coefficients for 
various frequencies, as needed for 
instance for flutter analysis. An 
example of the time saving obtained 
by using the two above features is 
given in the following section. 

^ 6. Numerical Results 

Typical numerical results obtained with 
EOUSSA (Ref. 14) arc presented in this sec- 
tion. Figure 1 shows the sectional lift 
coefficients at various stations of a wing- 
body in steady, subsonic flow compared 
against experimental and theoretical 


1 results of r.et. C and 15. The result^ 
lw:ri- eb*. nined for M=0 and i\ rcctanqul^’.r 
Jwing with chord c"l, and span b-6 thick- 
ness tr^O.09 and , The body io at 

zero angle of attack with overall length 
of 5 chords (lorcbody with length 
fuselage lergth L|;-3) . Note that . . 

tlie fu.'^r-lago ?.s cloned at the «’nd by n 
circ-.'lar piiito. In addition, a flat 'wake 
is e.T.anating frv.m the wing trailing edge 
and a cylindrical wake from that of the 
fu."clace. Figure 2 shows the convergence 
analyris of sectional lift of Figure 1 as 
a function of the number of elements. The 
corrut'-ir tir.o used for each cnr.o is-also 
indicated. Figures 3 and 4 present tho 
lift and r.iorent coefficients of a rec- 
tangular v.'ir.g in supersonic unsteady 
flows with aspect rat io. AR“2, c<>=l, b»2, 
T«0.001 and complex reduced frequencies 
k --^-0. 2+il.O, O.O+il.O and +0.2+11.0 
^ V 

for Each .unbor H=2 •' to 2.5. Tho re- 
sults an compared with those of Ref. 16. 
Noted that in contrast to Ref. 6, the 
present ret hod docs not require the use 
of diaphray.ms. Figure 5 prosentr; a 
wir.g-’ cdy-tail conHguration in fully 
unsteady ■^iov/ with specifications of the 
gcc.T.etry similar to that of Pig. 1. 

However, a h'r:.izontal tail is added with 
, chord c»l and b»6 which is statio.ned 0.5 
chords above the center line of the 
: fuselage. The complex reduced frequency 
I was k *C.l+i0.5. No existing result is 
available fer comparison, for, nr. men- , 

, tioned above, the present method is the 
only existing one which can analyze fully 
unsteady flow. The result is presented 
to demonstrate the generality of the 
method and its ability to handle fully 
unsteady flcv; problems. Figures C and 7 
presents the lift and moment coefficients 
and their corresponding phase angles for 
a rectangular wing oscillating in plunge 
and pitch in subsonic flov; with AR“2, 
T^O.COl, '.’ach number .’4x0. and 
4x7x7 elements on the whole wing. The 
results are identical with the ones 

i obtained \;ith SOSSA. However, consider- 
j. able time saying was obtained in the 
respective frequency and mode calculations 
; by the decomposition of the matrix Q into 
j frequency and node dependent matrices Qt 

I Jsee Sg . (43 ) )^. All the results wor e ob- 

I tained in 44 mins. 'i.e. 82.5 sees, for one 
aerodyna.mic coefficient and one frequency 
(since four coefficients and eight fre- 
quencies have been considered) . Note 
that the time for one single coefficient • 
and one single frequency (if evaluated 
independently) is about 13 minutes. 

Tables 1 and 2 contains tho general- 
ized forces for an AGARD wing-tail con- 
figuration in quasi-stcady and oscilla- 
tory flcv; compared with several existing 
. netheds (Refs. 17 to 21). While Table 3 
included the generalized forces for the 
■ same c o.a figuration in fully unsteady 
I flow (complex f requcncV)T“ Por'iriT' tlie're^' 
I suits the standard AGARD geometry (dcs- 
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crlbed for instance in Rofn. 18 end 19) was 
uecd. T lia consists of two cvept tnporcd 
lilting bur faces, ^ho first surface has 

and x^^i. 2.25 at y«0 and Xjj,»2.75 and 

x^^*3.70 at y»l and is located at r«0. The 
second surface has Xj^^-2.70 and x^j^^,«4.00 
I at y*0 and Xj^g«3.D0 and ct y^l 

land is located at zf.ri. All the renultc 
proscribed hero v?ere obtained using 4x7x7 
olaments on each surface. Results ob- 
I tained with 4x5x5 elen.ents indicate that 
j convergence was attained. The results are 
I usually in excellent agreewent with those 
I of Refs. 17 to 21^ ' 

p - . - 

7. Conclusions 

A general f emulation and computer 
program for the analysis of steady, oscil- 
latory and unsteady, subsonic and super- 
sonic aerodynamic flows around coupJox 
co.nfigurations have been presented. The 
final output of the code is the matrix 
of the aerodynamic influence coeffi- 
cients for flutter analysis to be used 
for instance in U'.e program FCAP (Ref. 

22 ). 

It should be noted that, while there 
exists several methods to analyze the 
problem of unsteady conpressiblo flows 
for complex configurations, the present 
method, embedded in the computer program 
SOUSSA, is unique in the following 
aspects t 

1. It provides a completely unified 
approach foi steady, oscillatory 
and fully un. teady, subsonic and 
supersonic ae,. .'dynamic flows. 

2. It can be applied to 'rbitrarily- 
complex configurations. V7ing-bo ay- 
tail configurations in fully ur- 

. steady flows hove been presented. 

3. It is computationally extremely 
general, flexible, efficient and 
above all, accurate. The elimina- 
tion of diaphragms in supersonic 
flow improved considerably the 
simplicity and efficiency of the 
code. 

4 . SOUSSA is the only existing pro- 
gram that can analyze fully unsteady 
complex-configuration potential 
aerodynamics in subsonic or super- 
sonic regimes. It is also the only 
program capable of handling oscil- 
latory supersonic aerodynamics for 
complex configurations. 


6. Flutter and opl<m.il design ana lyse a 
' rcquiic ev»iluntion of the aerodynamic 
Influence coefficients for several 
frequencies and mode shapes. With 
the unique features mentioned above, 
(i.e., separation of the aerodynamic 
influence coefficient m.ntrices into 
frequency and node dependent and 
independent matrices) the cotnputc’r 
time* that norm.illy would have been 
required is drairatically reduced. 

Appe n dix A 

In this Appendix the for.mulation for 
the suyarsonic case is briefly outlined. 
For conciseness, only supersonic trailing . 
edges are connidored ro that the contribu-, 
; tion of the wa)<o can ba ignored. (In j 
, SOUSSA diaphagms are not used end thcre- 
■ fore Ui 3 supersonic Wv'Jtc is treated as : 
' the sU'^ronic one.) Under snvall-petturba- ■ 
I tion assumption, U.c Green th-.orcM for 
potential supersonic flow is given by • j 


1-H-dr 




(a;i) 


where 


a 

=- ( is the conormal 

derivative .(Ref . 4) ) is the conormal 
wash which is prescribed by the boundary 
conditions, and * 

.Y»y/f‘ Z-r/x (A. 2 ) 

with Furthermore, 

« • [pt - X.)* - (Y - Y.)’ - (Z - Z.)’ ] (A. 3 ) 

whereas .... . . 

H -1 for X. - X> j^(y - Y*)^+ (Z - Z.f 

r 9 - 11 / 2 (A. 4 ) 

- 0 for X, - Xil CY - Y.r + CZ - Z . 


and 


• ®* 

[]-I 


(A. 5) 


T-e" 


5. In contrast to existing methods, 
which in many instances requires 
extensive user's background in 
aerodynamics and familiarity with 
the specific method, the present 
code requires very limited hu.T.an 
intervention and is, extremely easy 
to use. 


• I 

0* - M (X« - R* 

Following the sane procedure used for 
the subsonic case one obtains 
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!#•» V» s* ‘ 4 


n 

n - • . 

♦S‘>;s[\tr-3,*)-\(r-0j-)] 


i;: 




Tht definition of and is 

discussad at tha and of this appendix. 

Finally, taking the Laplace transforin 
of equation (26) results in 


yield] no problem for the elements 
com^>letely outnid*. or completely incido 
thu lUich forecui.e. However for iho 
elcr.ents partially inside the tSach 4 . 
forocone, a special dofinition for 9|^- 

and ^ nunt be uiodi note that 

and represent the two propagatibn 

tines for the dlcturbancan emanating from 
the olei.-ent to reach the point . 

Therefore for* alenonta^partially inside 
tha Koch foreconc, 0 , and g.r are most 

In* * jn 

appropriately defined '■(fron a physical 
point of viow) 06 the propagation Jkines. 
from the centroid of the portion of th<^ 
elerr.ei.t intersected by the Mach 

forecon^'to the control point Pj. 

It should be noted that with this 
definition of ,Q.. and oT. , supersonic 
flov;s can be trertcc! exaefi'y the samo 
way as suhronic flo..'s. rot instance 
a wing \/ith supersonic leading edge is 
solved by using both sides of the wit^g_ • 
sinultatKouoly. Also for wings' vrith 
partially supersonic leading edges the 
use of diaphragms is not necessary. 

This io a cor.sifiernblc advantage: the use 
use of diaphragr.E in the program SOSSA 
was cumbersome, especially for wing- 
body-tail and non-coplanar surfaces 
analyses. In the program SO’JCSA there 
is no differe.ncc in the treatment of 
subsonic and superoonic flows. 


where 


(e * - 1 


(A.IO) 


(A.ll) 


Equation (A. 9) yields the matrix 
to bo used in Eq. (39) for supersonic"... 
flows. The supersonic matrices , 

and are equal to the subsonic 

ones. Therefore the above modifications 

in the definitions of Yi. and z'.. are 

jh 3h 

the only changes necessary for supersonic 
flows. As mentioned above if the trailing 

edges are partially subsonic thie wake iB~| 
treated as in the subsonic case. ! 


Finally a fow remarks are necessary 
on tho fact that diaphragms are not used 
in SOUSSA and on how this ,relates_to the 
correct definition of and 'e.r It 


correct definition of 07 and e„ 

+lh jT**'. jn • 
should be noted that 8 |^' and 6 ^^ 1 

imaginary, if the centroid of the 


element 


lies outside tho Mach fore- 
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rig. 3. Litt coafficient, C, « for 
rectangular wing oscillatin9 in pitch, 
with AR-2, t-0.001, K “8, K -7. Co.n- 
parison with reaulta Of Raf? 16. 
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Fig. 4. Konent. courficient, , versus 
M, for rectangular wing oscJlXating in 
pitch, for AR“2, T-0.001, W -8, N "7. 
Cor.parison with results of ^ Ref. ^16. 



Fig. 5a. Pressure distributions at 2y/b 
■0.78 wing spanwise station of a wing- 
body-tail configuration in fully un- 
steady pitching mode with pitch axis 
at wing aid-chord. Conplex froque.'.cy 
k -0.l4il.5,Bpan b-6, furelage radius 
r°-0.5, thickness ratio T-. 09., total 
number of elements NBLEK - 388. 
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rig. 5b. Prcsfur* distributions at 2y/b 
■0.78 horizontal-tail sp&nwisc sLulion 
of the same configuration as Fig. So. 
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Pig. 7. Moment coefficient, C.^ , versus k 
for a rectangular vin? oscillating in 
plunge, with AR*2, t*0.001, M>0, N *7, t: 
7, N *20, wake length L "2c. Com parison ^ 
with^results of R o£. r— 


Pig. 6. Lift coefficient, , versus k, 
for rectangular wing oscillating in plunge, 
with AR-2, T-O.OOl, K-0, N -7, N -7, H «20, 
wake length X^"2c. Comparison with results 
of Rof. 16. ^ 
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• 

•. * 


-O.COCS 

0.0917 

Ref. *21 


• 

• 

-0.057l‘ 

-0.0! 30 

Pro: ent 

Wing twist 

Wing bending' 

1.2 

-0.07<r.‘ 

0.0301' 

Ref. '20 

. 


-0.0234' 

C.OCOl 

Kol. 21 

• 


■ 


C.C5S0 

rsvs ,t 

wing ttr.ding 

'Wing I'criiirig 

2.2 

-o.r.7'.'/ 
■ O'.Olf / 
-9.0246 

9.'/-U7 

0.2464 

0.221B 

:-ef. 23 
Mcf. 21 
Pi c!.rnt 

Tail roll 

King henditig 

3.2 

-C.OCl 

-0.0715 

0.035C 

C.CC15' 

-0.BC12 

-0.0908 

Re£. 20 
Kef. 21 
Prcsn.it 

Tail- pitch 

'.ring bendl ;g 

<,2 

-O.OifC 

-0.0CC2 

0-.033E 

0.0485 
f .0104 
-0.0606 

Jn-f. 20 
r.c£. 21 
Present 

Tail roll 

Tail roll 

.3', '3 

c.o:o3 

0.2622 

Re£. 20 


• 

, 0 , O'/CO 
O.OiCJ 

0.3170 
. 0.2838 

F.oi. 21 
Present 

Tall pitch 

Tail roll ■ 

4.3 

0.0072 

0.C3C£ 

O.02C3 

0.UC-: 

V.2208 

0.2240 

Ref. 20* 
lie I. 21 
Present 

Tail roll 

Tail pitch 

■ 3,4 

014517 

0.4<30 

0.5000, 

C.1632 
0.23 08 
0,1874 

ReC. 20 
.nc£. 21 
Pt'er.iiit 

Tail pitch 

Tail pitch . 

4,4 

0.2iC5 
■ 0..31C^ 

0.2588 

.0.3010 

P.of. 20 
Ref. 23 

> • ‘ 


. 

0.3724 

0.2354 

• Presint 


TABLE 3a 

Generalized Acrodynanic Force CrcfCicients 
for AGAPO Ving-7«il Iritcrferc.iCe H«=0.8, 4z/L»0.0 



CcncrallzcJ Caused l>y 
Force la pressure ia 


Wing twist 
Wing bending 
Tail r611 
Tail. pitch 



TADLE 3b 

Generalized Aerodynamic Force Cotrfficinnts 
for ACAfO Wing-Tail Inter Jcrt'ncc M'3.0, Az/L*0.6 



Cenerallxed 
Force In 

Caused by 
pressure lo 

i.J 

Wing twist 

Wing twist 

1-^3, 

t Tall roll 

s Tail roll 

1+3 

Wing bending 

King twist 

2+4, 

t Tail pitch 

( Tail toll 

14-3' 


mm 

*^3 


0.0130 

D‘.4004 

0.0464 

0.2654 

0.2860 

0.2872 

















